209 


W 


NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California 


THESIS 


i5.i! 


INBAND  SCATTERING  FROM  ARRAYS  WITH 
SERIES  FEED  NETWORKS 


Seunghoon  Lee 
September  1994 

Thesis  Advisor:  _ David  C.  Jenn 

improved  for  public  release;  distiibutioD  is  unlimited. 


DTIC  ^yj.'^LITY  ny 


^4  J1  15  003 


REPORT  DOCUMENTATION  PAGE 


Form  Approved  0MB  No, 
0704 


Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  re^nse,  including  the  time  for 
renewing  instruction,  searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and  completing  and 
reiicwing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of  this 
collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington  headquarters  Services.  Directmate 
For  Information  Operations  and  Reports,  121S  Jefferson  Davis  Highway.  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the 
Office  of  Management  and  Budget,  Paperwork  Reductioo  Project  (0704-0188)  Washingttm  DC  20S03. 


.  AGENCY  USE  ONLY  (laewMMAl 

2.  REPORT  DATE 

September,  1994 

y  REPORT  TYPE  AND  DATES  COVERED 
Master’s  Thesis,  Final 


TITLE  AND  Sl^TITLE 

INBAND  SCATTERING  FROM  ARRAYS  WITH  SERIES  FEED  NETWORKS 


AITHOR(S) 

Seunghoon  Lee 


PERFORMING  ORGANIZATION  NAME(S)  AND  AODRESSaS) 
Naval  Postgraduate  School 
Monterey  CA  93943-5000 


SPONSORING/MONnORlNC  AGENCY  NAMEtS)  AND  ADORESS(ES) 


S.  FVNDINC  NUMBERS 


8.  PERFORMING  ORGANIZATION 
REPtMtT  NUMBER 


la.  SPONSORING/MONTTORING 
AGENCY  REPORT  NUMBER 


II  SITPLEMENTARY  NOTES  The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  official 
policy  or  position  of  the  Department  of  Defense  or  the  L'  S  Government 


12*.  DISTRIBVTlON/AVAlLABILm  STATEMENT  Approved  for  pubbc  release, 
'  stnbuUon  unlimited 


12b  DISTRIBUTION  CODE 
A 


13.  ABSTRACT 

Approximate  formulas  for  the  inband  RCS  of  an  array  with  series  feed  have  been  derived  The  formulas  are  based  on 
the  hypothesis  that  an  inadent  wave  exates  forward  and  backward  traveling  vmrs  on  the  main  line  The  ^iproximate 
formulas  are  in  good  agreement  with  results  obtained  using  scanering  matrices,  thocby  verifying  the  asnunpuons  made 
in  the  approximate  soluuon  Spikes  in  the  RCS  pattern  have  been  identifled  with  specific  scattering  sources  in  the  array 
The  parameters  affecting  the  level  and  location  of  the  lobes  have  been  noted  The  main  adv  antage  of  the  approximate 
method  relauve  to  the  scattering  matrix  method  is  its  speed  The  scattering  matrix  method  requires  that  a  matrix 
equation  be  solved,  and  its  size  increases  with  the  number  of  array  elements 


14.  SUBJECT  TERMS 

Senes  feed  antenna,  RCS.  Antenna  RCS.  Scattering  Matrix 


IS.  NUMBER  OF  PAGES 

50 


17.  SECURITY 

CLASSIFICATION  OF 
REPORT 


18.  SECURITY 

CLASSIFICATION  OF 
THIS  PAGE 


19.  SECURITY 

CLASSIFICATION  OF 
ABSTRACT 


18.  PRICE  CODE 


28.  UMnATIONOF 
ABSTRACT 


Unclassified 


NSN  7540-0 1-280- 5  yai 


Unclassified 


Unclassified 


Standard  Fora  298  (Rev  2-89) 
lYcacnbed  bv  ANSI  Sid  239-18 


Approved  for  public  release,  distribution  is  unlimited 
INBAND  SCATTERING  FROM  ARRAYS  WITH  SERIES  FEED  NETWORKS 


by 

Seunghoon  Lee 
CPT ,  Korean  Army 

B  A  .  Korea  Military  Academy.  March  1987 


Submitted  in  partial  fulfillment 
of  the  requirements  for  the  degree  of 

MASTER  OF  SCIENCE  IN  SYSTEMS  ENGINEERING 

from  the 

NAVAL  POSTGRADUATE  SCHOOL 
September.  1994 


Author 


Electronic  Warfare  Academic  Group 


ii 


AF  RACJ 


Approximate  formulas  for  the  inband  RCS  of  an  array  with  series  feed 
have  been  (terived  The  formulas  are  based  (m  the  hypothesis  that  an  incident 
wave  excites  forward  and  backward  traveling  waves  on  the  main  line.  The 
approximate  fmmulas  are  in  good  agreement  with  results  obtained  using 
scattermg  matrices,  thereby  verifying  the  ass'  /^ons  made  in  the 
a;^oximate  solution  Spikes  in  the  RCS  pattern  ikiv'  i  een  identified  with 
specific  scattering  sources  in  the  array.  The  parameters  affecting  the  level  and 
location  of  the  lobes  have  been  noted  The  main  advantage  of  tLe  approximate 
method  relative  to  the  scattering  matrix  method  is  its  speed.  The  scattering 
matnx  method  requires  that  a  matrix  equation  be  solved,  and  its  size  increases 
with  the  number  of  array  elements. 


Aeossslon  For  j 

■TIS  9RA&I 

DTIC  TAB 

□ 

Unannounced 

□ 

Just,  Mention _ 

By - 

Distribution/  _ 

Avfll lability  Cedes 


i - 

Av.r,  11  and/or 

*.*•'  1 

Spec lal 

li''' 

_ — li; 

TABLE  OF  CONTENTS 


1.  INTRODUCTION .  1 

n.  ANTENNA  AND  FEED  SCATTERING  ANALYSIS .  4 

A.  DEFINITION  OF  RCS .  4 

B  ANTENNA  SCATTERING .  6 

C  FEED  REFLECTIONS .  8 

D  DERIVATION  OF  THE  SERIES  FEED  RCS  FORMULAS .  13 

III.  COMPUTATION  AND  RESULTS .  19 

IV.  CONCLUSION .  27 

APPENDIX  A.  COUPLING  COEFFICIENTS .  29 

APPENDIX  B.  MATLAB  PROGRAM  LISTING .  32 

APPENDIX  C.  SCATTERING  MATRIX  FORMULATION .  39 

LIST  OF  REFERENCES .  44 

INITIAL  DISTRIBUTION  LIST .  45 


IV 


I.  INTRODUCTION 


Current  state-of-the-art  radars  are  sensitive  enough  to  track  small  arms 
fire  and  even  a  swarm  of  insects.  In  response  to  this  sensitivity,  a  higher 
standard  of  stealthy  design  has  been  employed  to  avoid  detection  by  radar. 
The  result  has  been  radically  stealthy  platforms  like  the  F-117A  fighter  and 
B-2  bomber.  These  are  designed  to  have  the  smallest  possible  radar  cross 
section  (RCS),  low  infi'ared  (IR)  emissions,  and  low  microwave  emissions  in 
order  to  avoid  being  detected  and  tracked. 

The  RCS  of  a  target  not  only  de.oends  on  the  physical  shape  and 
material,  but  also  on  its  sub-components  such  as  antennas  and  other  sensors. 
These  components  on  the  platforms  must  also  be  designed  to  meet  low  RCS 
requirements  in  addition  to  their  sensor  system  requirements.  In  some  cases 
the  onboard  sensors  can  be  the  predominant  factor  in  determining  a 
platform’s  total  RCS.  A  typical  example  is  a  high  gain  anterma  on  a  low  RCS 
platform.  If  the  antenna  beam  is  pointed  toward  the  threat  radar  and  the  threat 
fi-equency  is  in  the  anteima  operating  band,  antenna  scattering  can  be 
significant. 

Scattering  fi’om  antennas  has  been  the  subject  of  interest  since  the 
1950s.  Extensive  work  has  been  done  with  regard  to  the  determination  of  the 
antenna  parameters  [Refl],  and  dipole  scattering  and  the  effect  of  the 
terminal  load  impedance  [Ref2].  The  RCS  of  boms  [Ref3],  reflector 
anteimas  [Ref  4],  and  microstrip  elements  and  arrays  [Ref  5]  have  also  been 
extensively  studied.  Detailed  analysis  of  the  inband  scattering  characteristics 
of  arrays  with  parallel  feed  networks  have  also  been  investigated  [Ref  6]. 
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This  thesis  examines  the  scattering  properties  of  phased  arrays  with 
series  feed  networks  (Figure  1)  which  operate  at  the  same  frequency  as  that 
of  the  illuminating  threat  radar.  For  this  so-called  inband  case,  the  threat 
signal  will  penetrate  into  the  feed  and  be  reflected  at  the  mismatches  and 
junctions  of  the  array.  This  occurs  even  for  a  well-matched  array  because 
small  mismatches  exist  within  the  feed  The  total  scattered  power  is  a  small 
fraction  of  the  incident  power,  but  for  a  large  array,  a  large  number  of 
reflections  exist  which  significantly  contribute  to  the  total  RCS  of  the  array. 
High  performance  arrays  must  not  only  meet  the  system  operating 
requirements  (i.e.,  high  gain,  low  sidelobe  levels,  etc.)  but  also  the  RCS 
requirements.  To  optimize  the  tradeoff  between  antenna  RCS  and  radiation 
performance,  an  efficient  and  accurate  model  of  RCS  is  crucial.  This  thesis 
presents  an  approximate  inband  scattering  model  for  arrays  with  series  feeds. 
The  approximate  model  is  shown  to  be  in  good  agreement  with  a  more 
rigorous  scattering  matrix  approach.  The  behavior  of  the  RCS  for  various 
feed  parameters  is  also  determined. 

Chapter  11  presents  the  theoretical  background  for  the  RCS  analysis  of 
the  series-fed  array  antenna,  and  contains  the  derivation  of  the  formulas  for 
the  approximate  solution.  Chapter  III  presents  some  computed  data  for  both 
the  approximate  and  rigorous  solutions.  Results  nic  also  presented  to  illustrate 
how  the  RCS  is  influenced  by  various  feed  parameters.  Finally,  Chapter  IV 
summarizes  the  results  of  the  data  acquired  and  concludes  with  some  design 
guidelines  for  controlling  the  behavior  of  the  RCS.  Some  conclusions  and 
recommendations  are  presented. 
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Pigarc  1 :  Series  feed  antenna  network. 


II.  A;4TENNA  AND  FEED  SCATTERING  ANALYSIS 


Scattered  fields  from  a  body  are  set  up  by  a  combination  of  induced 
electric  and  magnetic  currents.  The  currents  are  induced  such  that  the  fields 
satisfy  the  boundary  conditions  and  Maxwell’s  equations.  The  scattered  field 
can  be  determined  from  the  current  using  the  radiation  integrals.  However,  in 
most  practical  cases,  determining  die  induced  current  is  difficult  or 
impossible.  For  a  target  such  as  an  antenna  the  prediction  of  RCS  is  difficult 
because  of  many  scattering  sources  both  at  the  aperture  and  in  the  feed 
circuit.  In  fact,  even  antennas  with  identical  amplitude  and  phase  radiation 
patterns  can  differ  in  the  way  they  scatter. 

A.  DEFINITION  OF  RCS 

Radar  cross  section  is  a  measure  of  power  scattered  by  the  target 
towards  the  illuminator.  The  definition  of  RCS  can  be  stated  as 

power  reflected  to  receiver  per  unit  solid  angle 
incident  power  density  /  4n 

— ►  — ^ 

In  terms  of  the  incident  and  scattered  electric  field  intensities  and  Es  the 
RCS  is 


CT  =  lim  4nR^ 

R-xx> 
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(1) 
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where  R  is  the  distance  from  the  target  to  the  observation  point.  In  subsequent 
discussions  only  the  monostatic  RCS  is  considered.  That  is,  the  observation 
point  is  coincident  with  the  transmit  source.  The  limiting  process  in  (1) 
assures  that  the  incident  wave  is  approximately  a  plane  wave.  Calculation  of 
RCS  is  essentially  a  matter  of  finding  the  electric  field  scattered  from  a  target. 
If  the  current  induced  on  the  target  by  the  incident  plane  wave  can  be 
determined,  then  the  radiafion  integrals  used  in.  noatenna  analysis  can  be 
applied  to  compute  the  scattered  field.  However,  determining  the  induced 
current  is  a  difficult  problem  because  Maxwell’s  equations  must  be  solved  for 
complicated  boundary  conditions. 

The  unit  of  RCS  most  commonly  used  is  decibels  relative  to  a  square 
meter  (dBsm)  which  is  computed  by 

a,  dBsm  =10  log(a,m^).  (2) 

Radar  cross  section  depends  on  many  fectors  such  as  frequency  and 
polarization  of  the  incident  wave,  and  the  target  orientation  relative  to  the 
radar.  RCS  levels  can  range  from  10,000  m^  for  ships  to  less  than  0.001  m^ 
for  insects. 

The  scattering  characteristics  of  a  target  are  strongly  dependent  on  the 
frequency  of  the  incident  wave.  There  are  three  frequency  regions  in  which 
the  RCS  of  a  target  is  distinctly  different.  They  are  referred  to  as  the  low, 
resonance,  and  high-frequency  regions,  where  low  and  high  are  defined 
relative  to  the  size  of  the  target  in  terms  of  incident  wavelength,  rather  than 
their  physical  size.  For  a  target  of  characteristic  length  L,  the  low-frequency 
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or  Rayleigh  region  defined  by  kL  «  1  (  k  =  In/X).  The  second  is  the 
resonance  or  Mie  region  where  kL  »  1 .  The  last  is  the  high-fi-equency  region 
or  the  optical  region  defined  by  kL  »  1.  Depending  on  the  threat  radar’s 

wavelength  relative  to  the  antenna’s  operating  band,  an  arbitrary  target  can 
fall  into  either  of  the  three  regions. 


B.  ANTENNA  SCATTERING 

The  antenna’s  impact  on  RCS  is  twofold.  First,  die  antenna  fi’equently 
disturbs  the  continuity  of  the  surface.  This  introduces  more  edges  with  a 
potential  for  increasing  wide-angle  scattering.  A  second  aspect  of  the  antenna 
scattering  is  that  the  threat  radar’s  signal  can  penetrate  into  the  feed  and  be 
reflected  at  internal  mismatches  and  junctions.  Even  though  these  mismatches 
are  small  in  the  operating  band,  there  can  be  a  large  number  of  scattering 
sources  that  add  constructively  under  some  conditions.  This  scattering  is 
dependent  on  the  type  of  feed  and  the  devices  incorporated  therein. 

The  basic  equation  of  antenna  scattering  has  been  presented  by  several 
authors  [Ref  2,  7-9].  The  total  scattered  field  for  a  linearly  polarized  antenna 
when  the  port  is  terminated  with  a  load  is 


EsiZi)  =  E^iTa)  + 


/n 


— >  — ^ 
h{h 


(3) 


This  discussion  parallels  that  given  in  [Ref  6]. 
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where 


Z,  =  R,  +  jX,  is  the  radiation  impedance 

R,  =  +  Rj  the  antenna  resistance 

Rj  =  the  radiation  resistance 

Rj  =  the  ohmic  resistance 
— ^ 

h  =  effective  height 
— > 

Ei  =  incident  field 

T|  =  ^ji/e  is  the  intrinsic  impedance  of  the  medinm  and, 
r 

^  Zl  +  ZI  • 


(4) 


When  the  load  impedance  is  the  complex  conjugate  of  the  radiation 

impedance,  is  called  a  conjugate-matched  load.  In  equation  (3)  £t(Z*)  is 
the  structural  mode  and  the  second  term  on  the  h^-4iand  side  is  the  antenna 

mode.  When  the  antenna  is  conjugate  matched,  fise  antenna  mode  vanishes 
because  the  reflection  coefficient  T  goes  to  zero.  This  can  usually  be  achieved 
by  having  real  values  for  the  antenna  impedance  (  i.e.,  “tuning  out”  its 
reactance),  and  connecting  the  tenninals  to  a  matched  transmitter  or  receiver, 
thereby  having  no  reflection  at  the  junction  between  the  antenna  and  the 
transmission  line. 

Kahn  and  Kurss  [Ref  10]  found  that  for  a  dipole,  conjugate  matching 
results  in  half  of  the  power  being  scattered  and  half  absorbed.  The  reason 
being  that  for  an  isolated  dipole,  the  structural  mode  is  significant.  If  low 
RCS  is  desired,  then  it  is  necessary  to  mismatch  the  dipole  terminals  to 
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generate  sufficient  antenna  mode  to  cancel  the  structural  mode.  Forcing  the 
antenna  mode  term  to  zero  by  forcing  to  zero  is  impractical  for  large  arrays 
because  Z,  depends  on  angle  and  varies  with  element  location  in  the  array. 
However,  it  is  possible  in  principal  to  reduce  an  array’s  RCS  to  low  values. 
Green  [Ref  11]  has  shown  that  if  a  lossless,  linearly  polarized  antenna 
absorbs  more  power  than  it  scatters  under  matched  conditions,  its  gain  in  the 
backward  direction  must  exceed  the  forward  direction  gain.  This  condition  is 
satisfied  in  the  operating  band  of  high-gain  arrays. 

C.  FEED  REFLECTIONS 

In  the  case  of  an  operational  antenna,  the  incident  threat  signal  power 
delivered  to  the  radiating  element  terminals  actually  enters  the  feed  network. 
The  signal  transmitted  through  the  aperture  elements  is  reflected  at  points 
inside  of  the  feed  circuit.  For  a  phased  array,  the  feed  network  can  be 
extremely  comphcated,  with  potentially  every  transmission  line  discontinuity 
a  source  of  reflected  signals.  In  this  thesis,  the  signal  behavior  in  a  series  feed 
is  examined. 

A  series  feed  is  one  in  which  the  power  to  each  element  is  tapped  off 
of  a  main  line  sequentially,  as  fllustrated  in  Figure  1 .  The  coupling  values  are 
adjusted  to  provide  the  desired  amplitude  distribution.  The  most  common 
method  of  terminating  the  main  line  is  matched  loading;  another  is  to  use  a 
short  which  is  called  a  standing  wave  feed.  For  the  loaded  case,  some  power 
is  always  reflected  fi’om  the  load  at  the  end  of  the  main  line  which  gives  rise 
to  a  “reflection  lobe”  in  the  radiation  pattern.  It  is  reduced  in  magnitude 
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relative  to  the  main  beam  by  an  amount  determined  by  the  load  reflection 
coeflScient. 

There  are  two  antenna  scattering  modes  to  consider  in  the  calculation 
of  RCS.  The  first  is  the  antenna  or  radiation  mode  which  is  determined  by  the 
radiation  properties  of  the  antenna  and  it  vanishes  when  the  anteima  is 
conjugate  matched  to  its  radiation  impedance.  The  second  mode  is  the 
structural  mode  that  is  generated  by  the  induced  cmrents  on  the  antenna 
surfaces.  This  thesis  only  examines  the  antenna  mode  which  is  the  dominant 
RCS  contribution  to  a  phased  array  in  the  operating  band. 

Total  RCS  can  be  computed  by  summing  the  scattering  fi'om  all  the 
signals  that  enter  the  array,  are  reflected  by  the  mismatches  of  the 
components,  and  then  return  to  the  aperture.  The  RCS  is  obtained  fi’om 
combining  equations  (1)  and  (3)  over  all  the  array  elements  [Ref  6] 


a(e,<|)) 


4TlAe 

~ 


2 


(5) 


where 


r„  (0,(t))  =  total  reflected  signal  returned  to  the  aperture  for  element  n 

when  the  wave  is  incident  fi’om  the  (0,  <}>)  direction 

^  =  it(ux  +  \y  +  M'z)=  wave  vector 
w  =  sin  0  cos  ^ 

V  =  sin  0  sin  (|> 
w  =  cos  0 

d„  =  position  vector  to  element  n 

=  effective  area  of  an  element  =  h^Zgl^Ra 

^norm  ~  Dormalized  element  scattering  pattern. 
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There  are  two  approaches  to  obtaining  T  in  equation  (5).  One  is  to  use 
a  network  matrix  method  such  as  scattering  parameters.  All  interactions 
between  the  feed  devices  are  included.  The  scattering  matrix  iq)fMroach  is 
computationally  intense  because  a  matrix  equation  must  be  solved.  As  N 
increases  the  size  of  the  matrix  increases.  A  second  approach  is  to  trace 
signals  through  the  iced  However,  tracing  the  signal  flow  inside  a  feed 
network  can  be  extremely  difiicuh  and  tedious  job,  and  to  include  all  the 
multiple  reflections  for  accurate  calculation  is  even  more  difficult.  For  inband 
cases  where  the  reflection  coefficients  of  the  devices  are  small  (r  «  1), 
higher  order  reflections  can  be  neglected  The  following  approximations  are 
made  for  the  simplified  analysis  [Ref.  12] : 

1 .  Devices  of  the  same  type  are  assumed  to  have  identical  electrical 
characteristics.  For  example,  all  of  the  radiating  elements  have  the 
same  reflection  and  transmission  coefficients  (  and  ),  with  the 
exception  of  phase  shifters  which  have  a  transmission  coefficient  of  the 
form 


hn  = 


where 


Xn  =  («-l)a„ 

a,  =  ^o<f;tsin05cos((»5  =  kodxUs, 

(0,»  <l>s)  =  tfie  direction  of  the  array  radiation  beam. 
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Note  that  the  magnitude  of  the  transmission  coefficient  is  the  same  for 
all  phase  shifters  ;  \tp,  |  =  fp 

2.  The  feed  devices  are  well  matched  so  that  r  «  1  and  therefore 
hi^er  order  reflections  can  be  neglected. 

3.  Lossless  devices  will  be  assumed. 

With  the  previous  assumptions,  the  fraction  of  signal  entering  the 
radiating  elements  that  is  reflected  back  to  die  aperture  for  reradiation  is 


(6) 


where 


a  =  kji^u  is  the  incident  wave  inter-element  space  delay 
u  =  sin6cos(t>, 

E'„  represents  the  signal  re-emerging  from  the  series  feed  at  element  n.  As 
will  be  seen  later,  some  of  this  signal  is  due  to  elements  toward  the  input  (m  < 

n\  some  due  to  elements  toward  the  termination  (m  >  n),  and  some  due  to  the 
load  at  the  coupler  n. 
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Equati(Mi  (6)  shows  that  ^  tcrtal  severed  field  is  a  sum  of  scattering 
patterns  fi^om  each  mismatch  group  in  the  feed 

|£>|=l£;  +  £J+£J  +  £;|  (7) 

£5  is  the  total  scattered  field  due  to  the  series  feed  reflections.  It  is  obtained 
by  summing  over  all  as  described  in  the  next  section.  Assuming  that  only 
one  scattering  peak  dominates  at  any  given  incident  angle,  it  may  be 

sufficient  to  approximate  the  magnitude  squared  of  the  sum  in  (7)  by  the  sum 
of  the  magnitudes  squared 


|£1^»|£J|^+|£*|'  +  |£J|'  +  |£5|^  (8) 

The  RCS  contribution  fi-om  the  radiating  element,  phase  shifters  and 
coupler  inputs  are  [Ref  13] 


.  47l4^cos^6  r  sin(JVd)~|^ 

^  ^  _A^sin(a)_ 


(9) 


Gp 


Gc 


47L4^cos^e  2  fA  r  5in(Afa)~|^ 

p  *r  _jVsin(a)_ 


47l4^cos^9  2  fA  fA  sin(A^)  ^ 

h  ip  lNsmiQ_ 


(10) 

(11) 


where  ^  =  a  +  a, . 
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D.  DERIVATION  OF  THE  SERIES  FEED  RCS  FORMULAS 

Figure  2  shows  the  signal  distribution  for  the  n*  four  port  coupler  with 
coupling  coefficient  c„  and  transmission  coefficient  t,.  To  begin  with,  the 


Figure  2:  Coupling  and  transmission  paths  for  a  four  port  coupler. 

couplers  are  assumed  lossless,  matched  and  to  have  perfect  isolation  between 
the  coupled  port  and  through  port.  The  scattering  matrix  is 

0  t„  jc„  0 

fn  0  0  —jCn  ^j2'v 

jc„  0  0  t„  ^  ^ 

0  -jc„  t„  0 

To  obtain  formulas  for  the  feed  reflections,  the  received  incident  threat  signal 
can  be  decomposed  into  forward  and  backward  traveling  waves  on  the  main 
line.  Figure  3  shows  that  a  signal  received  by  element  n  appears  as  an 
excitation  for  the  N  -  n  elements  toward  the  load  and  the  n  -  1  elements 
toward  the  input.  Other  scattering  terms  not  included  in  the  forward  or 
backward  beam  are  the  “self-reflected  beam”  due  to  the  reflection  from  the 
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load  terminating  element  n,  and  the  *'input  load  reflected  beam”  due  to  the 
reflection  from  the  input  load  or  receiver  at  the  beginning  of  the  feed.  Thus 
for  each  element  in  the  array  four  beams  are  reradiated,  and  their  beamwidths 
and  peak  levels  depend  on  the  element  location  in  the  array.  The  total 
scattered  field  is  obtained  by  the  superposition  of  the  forward,  backward, 
self-reflected,  and  input  load  reflected  beams  from  all  elements. 


input  port 
loaded 

< - ► 


Figure  3  :  Forward  and  backward  traveling  waves  in  a  series  feed  network. 

The  formulas  for  the  RCS  of  the  series  network  can  be  obtained  by 
tracing  signals  through  the  feed.  Let  T,^  be  the  reflection  coefficient  of  all  of 
the  loads.  For  an  incident  signal  at  element  n,  “forward  wave”  that  travels 
toward  element  N  (see  Figure  4),  radiates  reflected  signals  from  n+1  up  to 
element  N  and  can  be  expressed  as 


(2)  Tbe  expressiOQs  for  the  electric  field  intensity  (E)  have  the  leading  constants  and  free  space 
Green's  function  omitted  because  only  a  ratio  is  needed  in  Equation  (1). 


14 


Ef^  =  eA"-lXa,4a) 

+  ^'0-  t„HjCn*2  +  ... 


+  /«l/«2  •  -  tN-l  jCN  ^xr-ixa,™)  }  _  (13) 

where  Vo~kL.  The  electrical  line  length  between  couplers  is  denoted  as 
and  it  is  equal  to  k  times  the  physical  length  of  the  line.  The  “backward 
wave”  that  travels  toward  the  input  (see  Figure  5),  radiates  reflected  signals 
from  the  n-1  elements  downline  from  element  n 


Figure  4  :  Forward  wave  (traveling  toward  element  N). 


+jc„-\e^'^- ^Ti  ° (oa-Kx) j 
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The  self-reflection  is  due  to  the  load  at  coupler  n,  and  it  is  unique  in  its  form. 
Thus  it  can  be  separated  and  simplified  as 


Figure  5  :  Backward  wave  (traveling  toward  element  1  -  input). 

For  the  reflection  from  the  input  load ,  the  equation  is  given  by 

n  (16) 

The  total  forward,  self,  input  load,  and  backward  traveling  fields  from 
all  the  elements  caused  by  an  incident  wave  is  obtained  by  summing  over  all 
elements 
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^-1  ( 

2 

f  m-1 

(17) 

M-l  1 

V|=1  J 

N-\ 

r  n-l 

(n-\  \  ] 

Eb  =  -TL  Z 

E  c 

(18) 

h^2 

L»m»1 

Vi«m  /  J 

(19) 

E„i/^  f  Ti  tl 

■ 

(20) 

The  minus  sign  occurs  for  some  terms  because  the  reflected  beams  pass 
through  two  coupled  arms  yielding  a  f  factor.  The  total  scattered  field  due  to 
the  series  feed  beam  is 

|£5|'  =  \Ef  +  Et+  E,„  + 

’«\Ej\^  +  \Ei\^  +  \Ej^  +  \Esei/  (21) 

The  correct  normalization  can  be  obtained  by  comparing  the  total 
scattered  field  in  equation  (21)  to  the  maximum  that  would  be  obtained  by  an 
array  with  completely  reflecting  elements  (rN).  Thus 


The  coupling  coefBcients  (cj  are  determined  by  the  prescribed  radiation 
beam  sidelobe  level.  Appendix  A  describes  how  the  coefficients  are 
computed  from  a  given  amphtude  distribution. 
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III.  COMPUTATION  AND  RESULTS 


Id  this  chapter,  calculated  results  are  presented  for  series  fed  arrays 
using  equations  (17)  through  (22).  To  verify  the  approximate  method,  data  is 
compared  with  the  rigorous  method  using  scattering  matrices.  The  scattering 
matrix  solution  is  described  briefly  in  Appendix  C.  The  aperture  element  is 
represented  by  a  two-port  device  with  the  reflection  coefficient  allowing 
direct  comparison  between  both  approximate  and  scattering  matrix  solutions. 
This  provides  a  means  of  evaluating  the  validity  of  the  assumptions  made  in 
deriving  the  approximate  equations. 

For  the  series-fed  phased  array  with  series  feed  netwcHk  operatmg  at 
wavelength  X,  calculations  were  performed  for  the  parametn^  sets  listed  in 
Table  1. 


TABLE  1  :  SERIES-FED  PHASED  ARRAY  PARAMETERS. 


parameter 

range 

no.  of  elements,  N 

50 

Vo 

0, 7t/4,  7i/2,  n 

d 

0.4  X 

0.2 

1 

0.5  X 

uniform, 

cosine-squared  on  a  pedestal 

e. 

0^45‘’ 
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A  MATLAB  program  for  calculating  the  approximate  RCS  is  shown  in 
i^pendix  B.  Each  element  of  the  array  (Figure  1)  that  is  illuminated  by  the 
incident  wave  has  its  own  forward  wave,  backward,  self  and  iuput  reflected 
term  as  discussed  in  Chapter  H,  Section  D. 

Some  parameters  are  varied  to  study  the  behavior  of  RCS.  They 
include  the  scan  angle  (0^ ),  the  line  length  between  couplers  (  ),  and  the 

coupling  coefficients  (  c„ ).  Figures  6  and  7  show  the  difference  between  the 
broadside  and  scanned  RCS  when  \|/^  =  7c/4.  It  can  be  observed  that  the 
highest  spike  at  6  =  0“  on  both  grajrfis  is  due  to  specular  scattering  from  the 
radiating  elements.  The  spikes  at  about  >20®  in  Figure  6  and  at  20®  in  Figure  7 
are  due  to  the  reflection  from  the  input  load.  The  spike  at  -35®  in  Figure  7  is 
due  to  Bragg  diffraction,  and  the  one  at  45®  is  due  to  the  effects  of  forward 
and  backward  waves  on  the  main  line  scattering  retrodirectively.  A 
comparison  of  the  scanned  and  broadside  cases  shows  that  the  reflections 


Figure  6  :  RCS  of  series  fed  array  by  the  approximate  method. 
(0s  =  O®,  v.  =  «/4) 
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from  points  behind  the  phase  shifters  move  toward  or  away  from  the  specular 
beam  along  with  the  array’s  radiation  beam. 


Figure  7  :  RCS  of  series  fed  array  by  the  approximate  method. 
(0s  =  Ji/4,  =  ic/4) 


In  general,  the  individual  spikes  can  be  attributed  to  specific  scattering 
sources  in  the  antenna.  The  specular  reflection  (  due  to  r, )  has  its  peak  at  9  = 
0°.  There  are  no  Bragg  lobes  for  this  term  because  the  element  spacing  is 
sufficiently  small.  Reflections  from  the  input  of  the  phase  shifter  have  a 
similar  behavior  as  can  be  seen  from  equation  (10).  Reflected  signals  from  the 
coupler  inputs  will  pass  through  the  phase  shifter  and  therefore  scan  along 
with  the  radiation  beam.  (This  can  be  a  problem  for  an  array  on  a  stealthy 
platform  if  it  is  used  to  track  a  threat.)  When  the  array  is  scanned,  a  second 
lobe  due  to  Bragg  diffraction  moves  into  the  visible  region. 


21 


30 


Figure  8  :  RCS  of  series  fed  array  by  the  approximate  method. 
(0s  =  0®,  v,  =n/2) 


Figure  9  :  RCS  of  series  fed  array  by  the  approximate  method. 
(0s=it/4,  V,  =ii/2) 
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A  large  lobe  also  arises  due  to  the  reflection  from  the  input  load.  In  this 
case,  the  array  acts  as  a  receive  antenna  that  collects  the  incident  threat 
signal.  A  portion  of  this  signal  is  reflected  at  the  terminals  and  appears  just  as 
if  a  transmit  signal  has  been  injected  into  the  antenna.  Thus  the  RCS  pattern 
is  the  gain  pattern  squared.  It  can  be  seen  by  comparing  Figures  6,  7  and  9 
that  the  lobe  from  the  input  load  reflection  is  determined  by  \\i^  as  well  as  0, . 


Figure  10  :  RCS  of  a  series  fed  array  by  the  scattering  matrix  method. 

(0.  =  O*,  V.  =  Jt/4) 

The  RCS  patterns  obtained  using  the  scattering  matrix  method  are 
shown  in  Figures  10  through  13.  They  correspond  to  the  same  cases  shown  in 
Figures  6  through  9.  This  solution  can  be  considered  rigorous  in  the  sense  that 
all  interactions  and  multiple  reflections  are  included.  The  agreement  between 
the  approximate  and  the  rigorous  solution  is  excellent. 
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Theta(deg ) 


Figure  1 1  :  RCS  of  a  series  fed  array  by  the  scattering  matrix  method, 
(0,  =  7C/4,  V.  =  5C/4) 


Figure  12  :  RCS  of  a  series  fed  array  by  the  scattering  matrix  method. 


Figure  14  :  RCS  of  a  series  fed  array  with  tapered  amplitude  distributiou  by 
the  approximate  method  ( 6, « ic/4 ,  y,  =  ic/2  ). 
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Figure  15  :  RCS  of  a  series  fed  array  with  tapered  amplitude  distribution  by 
the  approximate  method  (  6,  =  0*,  \|f,  =  ji/2  ). 

The  effects  of  changing  the  coupling  distribution  is  illustrated  in 
Figures  14  and  15.  The  RCS  behavior  is  essentially  the  same  but  the  feed 
contributions  to  RCS  are  reduced  slightly. 
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IV.  CONCLUSION 


i^proximate  formulas  for  the  inband  RCS  of  an  array  with  series  feed 
have  been  derived.  The  formulas  are  based  on  the  hypothesis  that  an  incident 
wave  excites  forward  and  backward  traveling  waves  on  the  main  line.  The 
approximate  formulas  are  in  good  agreement  with  results  obtained  using 
scattering  matrices,  thereby  verifying  die  assumptions  made  in  the 
approximate  solution. 

Spikes  in  the  RCS  pattern  have  been  identified  with  specific  scattering 
sources  in  the  array.  The  parameters  afifectiiig  the  level  and  location  of  the 
lobes  have  been  noted.  By  properly  choosing  parameters  such  as  \j/„  and  T^,  it 
is  possible  to  control  the  spikes  to  some  extent.  The  data  shows  that  there  is 
no  preferred  value  of  \\f^  that  eliminates  the  input  load  reflection.  It  is 
common  practice  that  the  line  length  between  couplers  should  be  set  to  an 
odd  multiple  of  a  quarter  wavelength  to  cancel  multiple  reflections.  However, 
total  suppression  of  the  RCS  can  only  be  achieved  by  perfectly  matching  all 
of  the  feed  devices,  which  is  practically  impossible. 

The  main  advantage  of  the  approximate  method  relative  to  the 
scattering  matrix  method  is  its  speed.  The  scattering  matrix  method  requires 
that  a  matrix  equation  be  solved,  and  its  size  increases  with  the  number  of 
array  elements.  A  typical  approximate  calculation  takes  about  2  min.  on  a 
Sparc  10,  whereas  the  matrix  solution  runs  about  2  hours  (50  elements  and 
0.5°  increments  ).  Also,  the  approximate  formulas  are  easily  extended  to  an 
arbitrary  number  of  elements.  Since  both  methods  assume  lossless  feed 
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networks  and  identical  reflection  coefScients  for  all  similar  devices,  the  actual 
value  of  RCS  can  be  expected  to  be  slightly  lower  than  shown. 

For  future  research,  it  is  recommended  that  the  RCS  of  a 
two-dimensional  array  be  investigated  and  be  presented  in  contour  form  [Ref 
6].  The  RCS  of  hybrid  feed  networks  (e.g.,  series  fed  in  one  dimension  and 
parallel  fed  in  the  other)  should  also  be  examined. 
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APPENDIX  A  :  COUPLING  COEFFICIENTS 


This  appendix  derives  the  series  feed  coupling  coefficients  for  a  given 
amphtude  distribution  {o„}.  For  uniform  illumination  (13dB  sidelobe  level) 
=1,  n=l,2,  N  where  N  is  the  number  of  couplers.  Also  of  interest  is  a 
cosine-squared  on  a  pedestal  distribution  (32dB  sidelobe  level)  where 


Figure  A.l  :  Cosine-squared  distribution. 
a„=ao  +  (l-  Qo)  sin^(^  x  7t) 

and  Go  =  0.2,  as  illustrated  in  Figure  A.I.  Referring  to  Figure  A.2,  for  a 
lossless  feed 

Pin  =11  kal+Pi 

n=l 

or. 
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Pin-Pi^k'Lal 

n 


and 


r.  _  f,.-n  fgr  ■  l-fl 
V,.2  ^  V^2 


Ho 

n 


Note  that  Pl  is  the  power  into  the  load  and  is  the  total  input  power. 


2  2 

kfl, 


Figure  A.2  :  Power  into  the  lossless  feed. 


The  coupling  coeflScients  are  defined  by 

""  p 

« 


where 


Pm 


m-l  _  m-1  _ 

H  kal  =  1-2  kal 

q=\  ^  9=1  ^ 
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Therefore,  for  unit  power  into  the  load  (P^  =1),  the  coupling  coefficient 
becomes 


Pin-  Z 
q=\ 


ai 


m-l 

-lo 

q=l 


N  «-l 


al 


n=l  w=l 


The  transmission  coefficients  are  given  by 


m 
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APPENDIX  B  :  MATLAB  PROGRAM  LISTING 


This  appendix  presents  a  listing  of  the  MATLAB  Code  that  computes 
the  RCS  using  the  approximate  formulas.  The  version  hsted  is  for  uniform 
amphtude  distribution  and  no  scan. 

%%  INBAND  SCATTERING  FROM  ARRAYS  WITH  SERIES  FEED  NETWORKS 
%%  Approximate  Method  (Broadside  RCS  Uniform  Distribution). 

N=50;,  Gamma=0.2;,  Si=pi/4; 

Ao=l;,  Pl=0.05;,  ki=0.8*pi*sm(0); 

theta=-pi/2:pi/360:pi/2;,  thetal=theta.*l  80/pi; 

%  a)  Forward  Wave  ;  (Travels  toward  element  N) 

for  a=l;N; 

A(a)=Ao+((l-Ao)*(sin(((a-l)/(N-l)).*pi)).^2); 

end 

for  1=1  :N; 

Cm(l)=sqrt((A(max(l))^2)/(((l/(l-Pl))*(sum(A(l:N).^2)))-(sum(A(l:l-l).^2))) 

); 

end 
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for  i=l;N-l; 

tl  (i)=((sqrt(l  -(Cm(max(i)).''2)))  *exp(j*Si)); 
t2(iHsqrt(l  -(Cin(max(i))."'2))); 


fork=l:N; 

tsf(k)=(l  -(Cm(max(i)).^2)); 
t3(k)=prod(tl(l:k-l)); 
t4(k)=prod(tl  (1  :k-2)); 
t5(k)=prod(t2(l  ;k-l)); 


for  z=l:361; 

del  1  (z)=sin(theta(max(z))).*(0.8*pi); 

del=dell(max(z)); 

for  n=l;N-l; 

for  m=n+l:N; 

Ef5(m)=Cm(m).*exp(j*(m-l)*(ki+del))*t3(m); 

end 

Ef4=sum(Ef5); 

EO(n)=Cm(n).  *exp(j*(n- 1  )*(ki+deI))*Ef4; 
end 

Ef2=sum(Ef3); 

Efl  (z)=-Gamma*Ef2; 
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%  b)  Self  term  ; 


forp=l:N; 

Esf3(p)=Gamma*tsf(p)  ♦exp(j*2*(p-l)*(ki+del)); 
end 

Esf2=sum(Esfi); 

Esfl(z)=(Es£2); 

%  c)  Backward  Wave  :  (Travels  toward  element  1 -input) 


for  q=2:N-l; 

for  r=l:q-l; 

Eb5(r)=Cm(r).*exp(j*(r-l)*(ki+del))*t4(r); 

end 

Eb4=sum(Eb5); 

Eb3(q)=Cm(q).*exp(j*(q'l)*(ki+del))*Eb4; 

end 

Eb2=sum(Eb3); 

Eb  1  (z)=-Gamma*Eb2; 

%  d)  Input  load  reflected  wave 

for  s=l:N 

Em3(s)=Cm(s)  ♦9xp(j*(s-l)*(ki+del+Si))*t5(s); 
end 
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Ein2=sum(Ein3); 

Einl  (z)=-Gamma*(l  *Gamma'^2)^2*(Em2).'^2; 
end 

Ef(:)=10.*logl0(abs(Efl(:))); 

Estl[:)=l  O.*logl  0(abs(Esn  (:))); 

Eb(:)=10.*logl0(abs(Ebl(:))); 

Em(;)=l  O.*logl  0(abs(Eml  (;);); 

Es=(abs(Efl)).^2+(abs(Esfl))."'2+(abs(Ebl)).^2+(abs(Einl)).'^2; 

Enl=Es./(N'^2); 

for  s=l:361; 

if  sin(dell(max(s)))~0 
Rf(s)=l; 
else 

Rf(s)=(sin(N*del  1  (max(s)))./(N.*sm(dell  (max(s))))).^2; 
end 

if  sin(0.8*pi*(sin(theta(max(s)))+siD<^0)))=0 
Rfs(s>=l; 
else 

Rfs(s)=(sin(N*0 . 8*pi*(sin(theta(max(s)))+sin(0)))./(N.  ♦sinCO .  8*pi*(sm(theta( 
max(s)))+sm(0))))).^2; 

end 

end 
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fort=l:361; 

RCSl(tH400*pi*(cos(theta(max(t))))^2).*(l-Gamma^2r4.*Enl(t); 

RCSr(t)=(400*pi*(cos(theta(max(t))))^2).*Gamma'^2.*Rf(max(t)); 

RCSp(t)=(400*pi'''(cos(theta(max(t))))'^2).*Gaimna^2*(l-Gamma'^2)'^2.*Rf( 

max(t)); 

RCSc(t)=(400*pi*(cos(theta(max(t))))'^2).*Gamma^2*((l-Gamma'^2)^4.*Rfs 

(max(t))); 

end 

RCS2(;)=1 0.*logl  0((RCS  1  (:))+ 1  .e-1 0); 
RCSR(;)=10.*logl0((RCSr(;))+l.e-10); 
RCSP(:)=10.*Iogl0((RCSp(;))+l.e-10); 
RCSC(:)=10.*logl0((RCSc(:))+l.e-10); 

RCSt=RCS  1  +RCSr+RCSp+RCSc; 

RCS(:)=10.*logl0((RCSt(:))+l.e-10); 

figure 

plot(thetal,Ef)5  g^d,  axis([-90  90 -50  30]); 
title(Torward  wave  vs  Angle’) 
xlabel('Theta(de,g.)'),  ylabel(Ef(dB)') 
pause 
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figure 

plot(thetal,Esf),  grid,  axis([-90  90 -50  30]); 
title('Self-Reflected  wave  vs  Angle') 
xlabel(Theta(deg.y),  ylabel('Esf(dB)’) 
pause 

figure 

plot(thetal  ,Eb),  grid,  axis([-90  90  -50  30]); 
titleCBackward  wave  vs  Angle') 
xlabel('Theta(deg.)'),  ylabel('Eb(dB)') 
pause 

figure 

plot(thetal,Ein),  grid,  axis([-90  90 -50  30]); 
titleCLoad  reflected  wave  vs  Angle') 
xlabel('Theta(deg. )'),  ylabel(Ein(dB)') 
pause 

figure 

subplot(22 1 ),  pIot(thetal  ,RCS2) 

grid,  axis([-90  90  -50  30]);,  title('Ant.  RCS  vs  Angle*) 

xlabel('Theta(deg.)').  ylabel(TlCS(dB)') 

pause 
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subplot(222),  pIot(thetalJlCSR) 

grid  axis([-90  90  -50  30]);,  titleCReflected  RCS  vs  Angle') 

xlabel('Theta(deg.)'),  ylabel(RCS(dB)’) 

pause 


subplot(223),  plot(thetal,RCSP) 

grid  axis([-90  90  -50  30]);,  title(Phase  Shifter  RCS  vs  Angle*) 

xlabel('Theta(deg.)'),  ylabel(RCS(dBy) 

pause 


subplot(224),  plot(thetal,RCSC) 

grid  axis([-90  90  -50  30]);,  title('Coupler  RCS  vs  Angle') 

xlabel(Theta(deg.)'),  ylabelCRCSCdB)') 

pause 


figure 

plot(theta  1  ,RCS),  grid 

axis([-90  90  -50  30]);,  title('Total  RCS  vs  Angle’) 

xlabel('Theta(deg.)'),  ylabel('RCS(dB)') 

save  rcsbn2  thetal  RCS  Ef  Esf  Eb  Ein  RCS2  RCSR  RCSP  RCSC 


pause,  axis;,  end 
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APPENDIX  C  :  SCATTERING  MATRIX  FORMULATION 


The  series  feed  is  represented  by  an  interconnection  of  two  and  four 
port  devices  as  shown  in  Figure  C.l.  The  unknown  quantities  are  the  set  of 
A  scattering  equation  can  be  written  at  each  port.  Thus  for  an  N 
element  array,  there  are  %N  equations  that  must  be  solved  simultaneously. 

The  scattering  matrix  of  the  radiating  element  is 


'  dn  dn  ■' 

r,  /,  1 

_  d2\  d22  _ 

— 1 

k. 

'Ki 

_ 1 

For  the  phase  shifter 

rp  ‘ 

The  excitations  are  due  to  the  incident  plane  where  S,  =  . 

Starting  with  element  1,  the  two  scattering  equatioas  at  the  radiating  devices 
are 


Pn  = 


(P")n  ^>^12 
I  (Pn)2i  iPn)22 


a\=S\d\i-^a3d\2  (1) 

02  =  S \d2\  +  (2) 
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and  at  the  phase  shifter  ports 


as  =  a2(pi)ii  ■^■^5(Pi)i2  (3) 

04  =  a2(pi)2i  +<^5051)22  (4) 

Four  equations  can  be  written  at  coupler  1  (cj 

06  =  00  rici)ii+\\f  09(01)12-^04(01)1^  +  01  r(ci)i4  (5) 

Oi  =  06  r(ci)2i  +M/  09^1)22  +  04(01)23  +  a?  r(ci)24  (6) 

05  =06  r(ci)3i  +11/09(01)32+04(01)33  +  07  r(ci)34  (7) 

01=06  r(Ci)4i  +V|/ 09(01)42 +04(01)43 +07  r(0l)44  (^) 

It  has  been  assumed  that  all  loads  are  identical  (r=  r,=  ...  =  F). 

Similarly,  all  interelement  line  lengths  are  equal  (v|/j=  ~ 

The  coupling  coefficients  are  given  in  Appendix  A  and  the  scattering  matrix 

by  equation  (12). 

For  elements  2  through  N,  a  pattern  repeats.  At  the  radiating  element 

08n-6  =  Sn  dn+  Oin-4  d\2  (8(n-l)  +  1) 

08n-5  =  Sn  0^21+  08n-4  d22  (8(n-l)  +  2) 

and  at  the  phase  shifter 

08«-4  =  08n-5(p«)ii  +  08n-2(Pn)i2  (8(n-l)  +  3) 

08m-3  =  Oi„-5(Pn)2i  +  0%n-2(Pi^21  (8(n-l)  +  4) 
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For  coupler  2  through  N-1 


a%n-l  =  w  ai„~i(Cn)ii  +  V|/  a%n+liCn)i2  +  ^8«-3(^^«)l3 

+^78«-ir(c„)i4  (8(n-l)  +  5) 

ain  =  Vj/  ain-sic„)2i  +  \\f  as„+l{c„)22  +  fl8»-3(^^«)23 

+ain-iTic„)24  (8(n-l)  +  6) 

ain-2  =  M/  ain-i(Cn)3\  +  M/  a^n+l(Cn)32  +  ^Sn-siCrds^ 

+ai„-\T(c„)24  (8(n-l)  +  7) 

Oin-l  —  ^8«-8(^«)41  +  ^8n+l(c^n)42  ^8/i-3(^^n)43 

+a%n-\TiCn)44  (8(n-l)  +  8) 

Finally,  coupler  N  must  be  treated  separately  because  of  the  extra  load 


a%s-l  =  Vj/  ^8A^-8(^A^)ll  +  CliN  T{CN)\2-^^iN-3iCN)i^ 

^a^N-xnctdn  («N-3) 

a%i^  =  V  ar87V-8(CA^)21+  r(CN)22  +  0%n-2(cn)23 

+ciiff-\T(cN)  24  (8N-2) 

OiN-l  =  v  aiN-i(cN) 31+  Cl%N  TiCN)32  +  a^N-3{CN)33 

+CI%N-\^{cn)  34  (8N-1) 

=  V  ^3r8A^-8(CAr)41+  Cl%N  r(CA^)42  +  87^-3 (cjv) 43 

+a8iV-ir(CA^)44  (8N) 


The  above  8N  equations  can  be  cast  into  matrix  form 


[S\  =  [F\[a] 


where  the  vector  [S]  contains  the  excitations  {S„}  and  [a]  the  unknowns  {a„}. 
[F]  is  a  feed  scattering  matrix  obtained  from  the  coefficients  of  the  above 
equations.  They  can  be  represented  as  entries  in  the  following  table 


column  =  unknown  index 


The  matrix  equation  is  solved  to  obtain  and  the  reflection  coefficients  r„ 
in  equation  (5)  can  be  calculated  from  ai ,  ajo, flgw-e,  •••,  ci%n-6- 
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